Budd-Vannimenus theorem for superconductors 
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The Budd-Vannimenus theorem is modified to the surface of superconductors in the Meissner state. 
This identity hnks the surface value of the electrostatic potential to the free energy in the bulk which 
allows one to evaluate the observed potential without the explicit solution of the charge profile at 
the surface. 



The Budd-Vannimenus theorem [1] links the value of 
the electrostatic potential at the surface of a metal to 
the energy per electron in the bulk. This simple conse- 
quence of the Feynman-Hellmann theorem turned out to 
be very useful in studies of metal surfaces [2] as it enables 
one to circumvent demanding calculations of the charge 
distribution at the surface. We will use this approach to 
evaluate the surface potential of superconductors within 
the Ginzburg-Landau (GL) theory. 

Perhaps it is useful to start from a historical review. 
Within the pre-London approaches based on the ideal 
charged liquid [3] , the electrostatic potential was assumed 
to balance the Lorentz and the inertial forces acting on 
currents. Since the current j = env flows along the sur- 
face and its amplitude falls off exponentially from the 
surface into the bulk, the electrostatic Bernoulli-type 
potential, e^p — — ^mw^, falls off on the scale of the 
London penetration depth. From the Poisson equation, 
p — — eoV^(/?, one finds that this corresponds to a charge 
accumulated in the layer penetrated by the magnetic 
field. Charge neutrality is maintained by the opposite 
surface charge. This picture has been confirmed by the 
London theory [4] . 

The electrostatic potential equals the Bernoulli poten- 
tial only at zero temperature when all electrons are in 
the superconducting condensate. At finite temperatures, 
a part rin of electrons remains in the normal state, while 
the rest = n — n„ contributes to the supercurrent, 
j = eUsV. As a consequence, the electrostatic potential 
reduces to et/? — —^^'m-v'^ as shown by van Vijfeijken 
and Staas [5]. Jakeman and Pike [6] recovered this re- 
sult from the time-dependent GL theory. Moreover, their 
approach describes the surface charge as a space charge 
localized on the Thomas-Fermi screening length. 

The surface charge and the charge accumulated at the 
London penetration depth form a surface dipole due to 
which the electrostatic potential in the bulk differs from 
the potential out of the superconductor. This voltage 
difference, which like the resistance of the Hall effect is 
perpendicular to the current, is not accessible by con- 
tact measurements but it was observed by a contact-less 
capacitive pickup [7,8]. 

The first measurements [7,8] were not sufficiently ac- 



curate and created expectations that the surface volt- 
age includes traces of the pairing mechanism. Rick- 
ayzen [9] used the thermodynamic approach and showed 
that the potential indeed has a pairing contribution, 
eip = -T^\mv^ - 4^4J^imw2. The first term re- 
suits from the Lorentz and inertial forces, while the sec- 
ond one reflects the pairing. The pairing term dominates 
close to the critical temperature T^- would be exper- 
imentally accessible, one could deduce from it the density 
dependence of Tc. The highly accurate measurement by 
Morris and Brown [10], however, chilled all the expec- 
tations. The electrostatic potential revealed no pairing 
contribution. 

The disagreement between the theory and experiment 
remained unexplained for a long time and the question 
of the charge transfer in superconductors was left aside 
till the discovery of the high-Tj, materials. It was pre- 
dicted [11,12] that in layered materials the superconduct- 
ing transition induces a charge transfer from Cu02 planes 
to charge reservoirs. This transfer caused merely by the 
pairing mechanism has been experimentally confirmed by 
the positron annihilation [13], the x-ray absorption spec- 
troscopy [14], and the nuclear magnetic resonance [15]. 

Apparently, there are two groups of contradictory ex- 
perimental results. The pairing contribution is absent in 
the surface potential but it is observed in the bulk. As it 
was indicated recently [16], there is an additional surface 
dipole which cancels the pairing contribution to surface 
potential seen in capacitive measurements. 

The approach presented in Ref. [16] employs the Budd- 
Vannimenus theorem, however, it has the same limitation 
as Rickayzen's formula - the magnetic field is described 
within the London's theory and it has to be weak not to 
perturb the density of the condensate. In contrast, the 
measurement of Morris and Brown [10] explores the en- 
tire range of the magnetic fields from weak fields up to the 
critical value. They found that the potential equals the 
magnetic pressure = -^2^1^^'^' independent 
of the temperature and the used material enters merely 
via the charge density put of its lattice, for both type-I 
and II superconductors. 

Here we modify the Budd-Vannimenus theorem so that 
it is consistent with the GL theory. Its validity is not 
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restricted to weak magnetic fields. We will show that the 
surface potential is given by 
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where /ei is the surface value of the GL free energy (3). 
With the help of the modified Budd-Vannimenus theorem 
one can explain the experimental result of Morris and 
Brown under very general conditions. 

The recent implementation of the GL theory has 
changed the picture of the surface charge. It is not local- 
ized on the scale of the Thomas-Fermi screening length 
as claimed by Jakeman and Pike, but it varies on scales of 
the GL coherence length [17]. As pointed out in Rcf. [18], 
if the electrostatic potential is a local functional of the GL 
wave function, (p[^], the GL boundary condition, V^" = 
in the direction normal to the surface, implies the zero 
electric field at the surface, E = —Vcp = f^VV' = 0. 
Accordingly, no additional surface charge is needed to 
maintain the charge neutrality. 

Since the surface charge is absent, it seems that the 
surface dipole on the Thomas-Fermi screening length is 
absent, too. There has to be yet another surface dipole, 
however. The pairing correlation is weaker on the sur- 
face than in the bulk, what results in forces pulling the 
Cooper pairs inside. Such forces are always balanced by 
the electrostatic field. The understanding of this effect 
will require microscopic studies which arc not yet feasi- 
ble. We can merely speculate that the surface dipole is 
somehow linked to the gap profile at the surface [19], i.e., 
it is localized on the BCS coherence length. 




^0 L 5 ^ X 

FIG. 1. Schematic picture of the electrostatic potential at 
the surface. On the scale of the London penetration depth 
A the potential ip{x) is concave, which corresponds to charge 
accumulation. Below the GL coherence length ^, the charge 
is doplotod loading to a convex potential. Within the BCS 
coherence length ^o, the potential makes a sharp step due to 
the surface dipole. (po is the true surface potential, while ^5(0) 
is the extrapolated GL value. The dashed line indicates the 
virtual compression of the lattice. 

To introduce the surface dipole on an intuitive level we 
first assume that the system is close to the critical tem- 
perature. In this regime, the London penetration depth 
A and the GL coherence length ^ are much larger than 
the BCS coherence length . We can then define an in- 



termediate scale L such that <^ L S,, X, as sketched 
in Fig. 1. On the scale L, the GL wave function, the vec- 
tor potential, and the electrostatic potential change only 
negligibly, e.g., ip{x) « (p{x ^ 0) = (p{0) for < ^ < L. 
The extrapolation of the bulk potential towards the sur- 
face, <p(0), has to be distinguished from the true surface 
potential ifQ. The difference ifQ — (p{0) is caused by the 
surface dipole we aim to evaluate. 

Close to the critical temperature we can follow the idea 
of Budd and Vannimenus [1]. Let us assume a virtual 
compression of the crystal lattice such that the lattice 
charge density is removed from the surface layer of an 
infinitesimal width SL. The perturbation of the lattice 
charge density in the infinitesimal layer < x < SL is 
Spiat = — Plat- The compression leads to an increase of 
the charge density in the layer SL < x < L, where (5/7iat = 
P\a,tSL/L, in accordance with the charge conservation. 

Let us remind the basic idea of the Feynman-Hellmann 
theorem. The lattice charge enters the jcUium model 
of metals as an external parameter. If one changes 
this external parameter, the situation corresponds to 
doing work on the system, SW = S J dx5p\g,t^^ = 
S J dxSpii,t<fi, where / is the density of the free energy 
including the electrostatic interaction, and 5 is the sam- 
ple area. According to the Feynman-Hellmann theorem, 
the change of the electrostatic potential does not con- 
tribute to the work up to the first order in 5p\at ■ Now we 
can proceed with the algebra. Wo split the integral into 
three parts. Since SL is an infinitesimal displacement, 
the potential in the layer < x < SL can be replaced 
by the surface value (pQ. The surface region SL < x < 
gives a negligible contribution of the order oi ^q/L. In 
the remaining bulk region < x < L, the electrostatic 
potential is nearly constant and equals (^(0). The work 
thus reads SW = SSLp\a,t{(f{0) — tpo). 

In the same time, the work increases the free energy T 
of the system, SYJ = ST . We denote the spatial density 
of the electronic free energy by /ei. Lower case letters 
denote free energies per unit volume, i.e., their densities, 
while the calligraphic upper case letters denote the cor- 
responding total values. The total free energy changes 
^^^=-'-W^SSL 



(-/el + SSL. The first 

term results from the reduced volume, SL S{L — SL), 

and the second one from the corresponding increase of 
the electron density by n — *■ n(l -|- SL/L). We thus ob- 
tain a modification of the Budd-Vannimenus theorem, 
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This relation is our final result. It describes the step of 
the potential at the surface due to the surface dipole in 
terms of the free energy. 

Formula (2) differs from the original Budd-Vannimenus 
theorem in three points. First, in the original Budd- 
Vannimenus theorem, the surface potential is related to 
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the potential ipoo deep in the bulk. In (2), the extrapo- 
lation of the internal potential towards the surface, </'(0), 
appears instead. Second, in order to cover systems at 
finite temperature, wc use the free energy instead of the 
ground-state energy. Third, within the original Budd- 
Vannimenus approach, the electron charge density and 
the lattice charge density have to be equal because of 
the charge neutrality. In our approach, the density of 
electronic charge differs locally from the lattice charge 
density, en ^ —piat, due to the charge transfer on the 
scales ^ and A. 

Let us demonstrate how the relation (2) can be used 
within the GL theory. To this end we introduce the GL 
free energy 



/el = 
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where ip is the GL wave function, A is the vector potential 
and TO* = 2to and e* = 2e arc the mass and the charge of 
the Cooper pair. The GL parameters a and /3 depend on 
the temperature and the electron density n = nn + 2\^p\'^ . 
Finally, we add the electromagnetic energy 

/ = /ei + ^(Aat+en)-^i?2 + -^i32, (4) 

with the magnetic field B = W x A and the electric field 
E = -V(p. 

Variations of the free energy with respect to its inde- 
pendent variables A, Lp,tp,nn yield the equations of mo- 
tion in Lagrange's form 
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For h' = A the variational condition (5) yields the 
Ampere-Maxwell equation, ior ly = (p the Poisson equa- 
tion, for v — ijj the GL equation, and for i/ = n„ the 
condition of zero dissipation. 
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This condition allows one to evaluate the electrostatic po- 
tential in the bulk of the superconductor [20]. Of course, 
one can add any constant to the electrostatic potential. 
The actual choice (6) simplifies the relations we deal with 
below. Another convenient choice would be to subtract 
the value of the potential in the non-magnetic state, i.e., 
deep inside the superconductor, so that <p approaches 
zero in the bulk. 

Formula (6) does not cover the surface dipole on the 
scale ^0, therefore at the surface it provides the extrapo- 
lated bulk value (p{0). We can thus use (6) to rearrange 
the Budd-Vannimenus theorem (2) as 



Platy'O = /el + V(0) (Plat + en) . 
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Now all terms on the right hand side are explicit quan- 
tities which one obtains within the GL theory extended 
by the electrostatic interaction [20]. 

In customary GL treatments, the electrostatic poten- 
tial and the corresponding charge transfer are omitted. 
This approximation works very well for magnetic proper- 
ties since the relative charge deviation, (piat + en) /piat, 
is typically of the order of 10~* to 10"'* leading to com- 
parably small corrections in the GL equation. With the 
same accuracy one obtains the electronic free energy /d. 
Therefore it is possible to evaluate the surface potential 
using (1), which follows from (7) if terms proportional to 
{put + en) /piat are neglected. 

So far we have restricted ourselves to systems close to 
the critical temperature when the validity conditions of 
the GL theory are well satisfied. In many cases, the GL 
theory is used beyond the limits of its nominal applica- 
bility, however. In these cases, the GL coherence length 
^ and/or the London penetration depth A are compara- 
ble to or even shorter than the BCS coherence length 
^0- Although the above derivation does not apply to this 
case, it is possible to show that the Budd-Vannimenus 
theorem (7) provides physically consistent predictions in 
these cases, too. 

We evaluate now the voltage which develops across a 
superconducting slab with the magnetic field parallel to 
the slab. For the slab geometry, the GL equation has an 
integral of motion, see Bardeen [21]. This integral can be 
obtained quite generally by the Legendre transformation 
of the free energy, g = f — 'Y^^{df ld'Wu)Wv. Indeed, if 
the fields v obey the equations of motion (5), the gradi- 
ent V(? = vanishes, i.e. g = const. With the help of 
T.Adf/dVy)Vv = (?iV™*)|VV'|' - eo-E^ + BVmo, we 
can express the electronic free energy as 
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<^(piat + en). (8) 



At the surface, the GL boundary condition demands that 
Vip = what implies E = 0. The free energy on the 
surface thus reads 



B 

fei= 9 + <p(0)(piat + en). 
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From the Budd-Vannimenus theorem (7) then follows 
Piat'Po = (2/xo), and since g = const, we obtain the 

potential difference across the surface as the difference of 
left and right magnetic pressures 
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Exactly this value of the voltage has been observed by 
Morris and Brown [10]. 

As noticed already by Bok and Klein [7], there is a 
simple argument for the formula (10). The left hand side 
represents the electrostatic force (per unitary area) on 
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the lattice, Foist = dxEput = Piat(</?o " '^o')- The 
right hand side is the Lorentz force FLor = BJ with the 
mean magnetic field B = ^{Bl+Br) and the net current 

J = dx j given by Ampere's rule, — Br = ^lqJ. 
Since the electrostatic field provides the only mechanism 
by which the force is passed from the electrons to the 
lattice, the two forces have to be equal, Flof = -Feist, 
which yields (10). 

Turning the argument around, one can view the deriva- 
tion of the Lorentz force as an alternative proof of formula 
(7). For the slab geometry, this is more general than the 
proof based on the Budd-Vannimenus approach, as the 
approach via the Lorentz force does not rely on the spe- 
cific order of the characteristic lengths ^, A needed 
above to define the intermediate length L. Moreover, the 
Legendre transformation provides the integral of motion 
for any modification of the GL free energy which is of the 
form /ei[V', Vt/)]. 

On the other hand, within the Budd-Vannimenus ap- 
proach the surface dipole is treated as a property of the 
superconducting condensate, what encourages us to hope 
that formula (7) or its approximation (1) can be used to 
obtain the surface potential also for cases when the mag- 
netic field has a component perpendicular to the surface. 
In particular, we expect that it will be applicable also 
to the superconductors in the mixed state, especially to 
evaluate the electric field generated by vortices penetrat- 
ing the surface [22]. 

In conclusion, the Budd-Vannimenus theorem was 
modified to the surface of a superconductor. It allows 
one to evaluate the electrostatic potential of the surface 
from the free energy and the bulk electrostatic potential 
near the surface. Within approximation (1) one obtains 
the surface potential from the free energy without the 
actual knowledge of the bulk potential. 
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